QCD at Finite Isospin Density 
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QCD at finite isospin chemical potential /ij has no fermion sign problem and can be studied on the 
lattice. We solve this theory analytically in two limits: at low fii, where chiral perturbation theory 
is applicable, and at asymptotically high fii, where perturbative QCD works. At low isospin density 
the ground state is a pion condensate, whereas at high density it is a Fermi liquid with Cooper 
pairing. The pairs carry the same quantum numbers as the pion. This leads us to conjecture that 
the transition from hadron to quark matter is smooth, which passes several tests. Our results imply 
a nontrivial phase diagram in the space of temperature and chemical potentials of isospin and baryon 
number. 



Introduction. — Ample knowledge of QCD in the regime 
of finite temperature and baryon density is crucial for 
understanding a wide range of phenomena from heavy 
ion collisions to neutron stars and cosmology. First- 
principles lattice numerical Monte Carlo calculations 
provide a solid basis for our knowledge of the finite- 
temperature regime. However, the regime of finite baryon 
chemical potential hb is still inaccessible by Monte Carlo 
because present methods of evaluating the QCD partition 
function require taking a path integral with a measure 
which includes a complex fermion determinant. Ignoring 
the determinant (as in the popular quenched approxi- 
mation) leads to qualitatively wrong answers for finite 
fig jjj- Such a contrast to the case of hb = 0, where the 
quenched approximation proved useful, comes from the 
fact that the latter corresponds to an unphysical theory 
with pairs of quarks of opposite baryon charges (conju- 
gate quarks) [0. This is one of the main reasons why 
our understanding of QCD at finite baryon density is 
still rudimentary. Many interesting phenomena, such as 
color superconductivity and color-flavor locking H, oc- 
cur at finite baryon density, beyond the reach of current 
lattice techniques. 

To understand the regime of finite baryon density one 
would need to follow the transition from hadronic to 
quark degrees of freedom by increasing the density of 
a conserved charge (such as baryon number), i.e., with- 
out invoking the temperature. This is the motivation 
for us to turn to QCD at finite chemical potential /// of 
isospin (more precisely, of the third component, 1$). Na- 
ture provides us with nonzero \xi systems in the form of 
isospin-asymmetric matter. These always contain both 
isospin density and baryon density. In any realistic set- 
ting \x\ <C fJ-B- In this paper, however, we shall consider 
an idealization in which /i/ is nonzero while /is = 0. 
Such a system is unstable with respect to weak decays 
which do not conserve isospin. However, since we are in- 
terested in the dynamics of strong interaction alone, one 
can imagine that all relatively slow electroweak effects 



are turned off. Once this is done, we have a nontrivial 
regime which, as has been emphasized recently in is 
accessible by present lattice Monte Carlo methods, while 
being, as we shall see, analytically tractable in various in- 
teresting limits. As a result, the system we consider has 
a potential to improve substantially our understanding 
of cold dense QCD. This regime carries many attractive 
traits of two-color QCD |||| , but is realized in a physi- 
cally relevant theory — QCD with three colors. 

Positivity and QCD inequalities. — Since the Euclidean 
version of our theory has a real and positive fermion de- 
terminant, some rigorous results on the low-energy be- 
havior can be obtained from QCD inequalities HI. In 
vacuum QCD, the latter rely on the following property 
of the Euclidean Dirac operator T> = 7 • (d + iA) + m: 



75^75 = 



(1) 



which, in particular, implies positivity det T> > 0. For the 
correlator of a generic meson M — ipTyj, we can write, 
using (Fy) and the Schwartz inequality: 



(M(x)Mt(O))^ - -(TtS(x,0)TS{0,x)T) a = 
(TrS(x,0)Ti l5 S\x,0)i l5 T) A < {TrS{x,0)S^{x,0)) A , (2) 

where S = D^ 1 and T = 701^70. The inequality is satu- 
rated for mesons with T = ij^Ti, since T> commutes with 
isospin Ti , which means that the pseudoscalar correlators 
are larger, point-by-point, than all other 7=1 meson cor- 
relators Q . As a consequence, one obtains an important 
restriction on the pattern of the symmetry breaking: for 
example, it cannot be driven by a condensate of {ip^ip), 
which would give + Goldstones. 

At finite isospin density, fii 7^ 0, positivity still holds 
[pj and certain inequalities can be derived (in contrast 
with the case of fis 7^ when there is no positivity) . Now 
T> = 7 • (d + iA) + hfJ>i7o T 3 + m i an d Eq. (|l|) is not true 
anymore, since the operation on the right-hand side of (|l|) 
changes the relative sign of [ij. However, provided m u = 
md, interchanging up and down quarks compensates for 
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this sign change (the u and d quarks play the role of 
mutually conjugate quarks 0]), i.e, 

ri7 5 2>75Ti =P f . (3) 

Instead of isospin t% in (||) one can also use (but not 
T3). Equation (||), in place of the now invalid Eq. (@), 
ensures that detl? > 0. Repeating the derivation of the 
QCD inequalities, by using (|j|) we find that the lightest 
meson, or the condensate, must be in channels i^i^Ti^, 
i.e., a linear combination of tt~ ~ u'jsd and ir + ~ c?7su 
states. Indeed, as shown below, in the two analytically 
tractable regimes of small and large fii the lightest mode 
is a massless Goldstone which is a linear combination of 
M75<i and dj^u. 

Small isospin densities. — When /i/ is small compared 
to the chiral scale (taken here to be m p ), we can use 
chiral perturbation theory. For zero quark mass and 
zero [ii the pions are massless Goldstones of the spon- 
taneously broken SU(2)l x SU(2)a chiral symmetry. If 
the quarks have small equal masses, the symmetry is 
only SU(2)l+_r- The low-energy dynamics is governed 
by the familiar chiral Lagrangian for the pion field E € 
SU(2): C = i/^Tr^E^Et - 2TO 2 .ReE], which contains 
the pion decay constant f„ and vacuum pion mass as 
phenomenological parameters. The isospin chemical po- 
tential further breaks S\J(2)l+r down to XJ(1)l+r- Its 
effect can be included to leading order in \i\ without 
additional phenomenological parameters by promoting 
SU(2)lxSU(2)a; to a local gauge symmetry and view- 
ing /ij as the zeroth component of a gauge potential Q . 
Gauge invariance thus fixes the way /ii enters the chiral 
Lagrangian: 

£ cff = ^LTrV.EV^t - I^ReTrE, (4) 

where the covariant derivative is defined as 

V E = a E-^(T3E-Er 3 ). (5) 

By using (|]) , it is straightforward to determine vacuum 
alignment of E as a function of /i/ and the spectrum 
of excitations around the vacuum. We are interested in 
negative fij, which favors neutrons over protons, as in 
neutron stars. The results are very similar to two-color 
QCD at finite baryon density B: 

(i) For < m^, the system is in the same ground 
state as at /i/ = 0: E = 1. This is because the lowest 
lying pion state costs a positive energy m v — to excite, 
which is impossible at zero temperature. 

(ii) When \fii\ exceeds it is favorable to excite tt~ 
quanta, which form a Bose condensate. In the language of 
the effective theory, such a pion condensate is described 
by a tilt of the chiral condensate E, 

E = cos a + i (ti cos 4> + t 2 sin </>) sin a , 

cos a. = tti^/ . (6) 



The tilt angle a is determined by minimizing the vac- 
uum energy. The energy is degenerate with respect to 
the angle 0, corresponding to the spontaneous breaking 
of the U(l)i + i? symmetry generated by I 3 in the La- 
grangian (|4|). The ground state is a pion superfluid, with 
one massless Goldstone mode. Since we start from a the- 
ory with three pions, there are two massive modes which 
can be identified with ttq and a linear combination of 7r + 
and 7T~. At the condensation threshold, m no = and 
the mass of the other mode is 2m,,-, while for \fij\ 3> 
both masses approach \fij\. 

The isospin density is found by differentiating the 
ground state energy with respect to /ij and is equal to: 

ni = fl^i sin 2 a = fim ^1 - , > m„ . (7) 

For just above the condensation threshold, — 
to^ <ti TOjt, Eq. (0) reproduces the equation of state of 
the dilute nonrelativistic pion gas ||. 

It is also possible to find baryon masses, i.e., the en- 
ergy cost of introducing a single baryon into the system. 
The most interesting baryons are those with lowest en- 
ergy and highest isospin, i.e. neutron n and A - isobar. 
There are two effects of /ij on the baryon masses. The 
first comes from the isospin of the baryons, which effec- 
tively reduces the neutron mass by and the A~ 
mass by Alone, this effect would lead to the for- 
mation of baryon/antibaryon Fermi surfaces, manifested 
in nonvanishing zero-temperature baryon susceptibility 
Xb = dns/d^B when fj,] > |mA- However, long before 
that, another effect turns on: the 7r _ 's in the condensate 
tend to repel the baryons, lifting up their masses. These 
effects can be treated in the framework of the baryon 
chiral perturbation theory Q , giving 

m n = rriN — cos a, m^- — m A ^ — cos a (8) 

in the approximation of nonrelativistic baryons. Equa- 
tion (||) can be interpreted as follows: as a result of the 
rotation (||) of the chiral condensate, the nucleon mass 
eigenstate becomes a superposition of vacuum n and p 
states. The expectation value of the isospin in this state 
is proportional to cos a appearing in ^ . With cos a 
given in Eq.([|), we see that the two mentioned effects 
cancel each other when m, <ti m p . Thus the 

baryon mass never drops to zero, and xb — at zero 
temperature in the region of applicability of the chiral 
Lagrangian. 

As one forces more pions into the condensate, the pions 
are packed closer and their interaction becomes stronger. 
When [ii ~ TO p , the chiral perturbation theory breaks 
down. To find the equation of state in this regime, full 
QCD has to be employed. As we have seen, this can be 
done using present lattice techniques since the fermion 
sign problem is not present at finite similar to the 
two-color QCD §. 
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Asymptotically high isospin densities. — In the opposite 
limit of very large isospin densities, or S> m p , the de- 
scription in terms of quark degrees of freedom applies 
since the latter are weakly interacting due to asymp- 
totic freedom. In our case of large negative /i/, or n.j, 
the ground state consists of d quarks and u antiquarks 
which, neglecting the interaction, fill two Fermi spheres 
with equal radii |/i/|/2. Turning on the interaction be- 
tween the fermions leads to the instability with respect 
to the formation and condensation of Cooper pairs, simi- 
lar, to some extent, to the diquark pairing at high baryon 
density 0. In our case, fii < 0, the Cooper pair consists 
of a u and a d in the color singlet channel. The order 
parameter has the same quantum numbers as the pion 
condensate at lower densities, 

(ufd) + 0. (9) 

Because of Cooper pairing, the fermion spectrum ac- 
quires a gap A at the Fermi surface, where 

A = 6|/i/|£T 5 e~ c/ff , c = 3vr 2 /2 (10) 

where g should be evaluated at the scale \ni\. This be- 
havior comes from the long-range magnetic interaction, 
as in the superconducting gap at large /is pj. The 
constant c is smaller by a factor of y/2 compared to the 
latter case due to the stronger one-gluon attraction in 
the singlet qq channel compared to the 3 diquark chan- 
nel. Consequently, the gap ( |l0| ) is exponentially larger 
than the diquark gap at comparable baryon chemical po- 
tentials. By using the methods of JPI| , one can estimate 
fewlO 4 . 

The perturbative one-gluon exchange responsible for 
pairing at large /i/ does not distinguish ud and S75C? 
channels: the attraction is the same in both. The u^d 
channel is favored by the instanton-induced interactions, 
which explains the fact that the condensate is a pseu- 
doscalar and breaks parity. This is consistent with our 
observation that QCD inequalities also constrain the 
1=1 condensate to be a pseudoscalar at any /i/. 

Quark-hadron continuity. — Since the order parameter 
(^) has the same quantum numbers and breaks the same 
symmetry as the pion condensate in the low-density 
regime, it is plausible that there is no phase transition 
along the ^/ axis. In this case the Bose condensate of 
weakly interacting pions smoothly transforms into the 
superfluid state of ud Cooper pairs. The situation is very 
similar to that of strongly coupled superconductors with 
a "pseudogap" [Q2J , and possibly of high-temperature su- 
perconductors |13| . This also parallels the continuity be- 
tween nuclear and quark matter in three-flavor QCD as 
conjectured by Schafer and Wilczek H]. We hence con- 
jecture that, in two-flavor QCD, one can move continu- 
ously from the hadron phase to the quark phase with- 
out encountering a phase transition. Since a first order 
deconfinement phase transition at intermediate isospin 



chemical potential cannot be rigorously ruled out (though 
it is unlikely, see below), this conjecture needs to be ver- 
ified by lattice calculations. 

A number of nontrivial arguments support the conti- 
nuity hypothesis. One notices that all fermions have a 
gap at large \m\, which implies that xb = at T = 0. 
This is also true at small /ij. It is thus natural to expect 
that xb remains zero at T = for all which also 
suggests one way to check the continuity on the lattice. 

Another argument comes from considering the limit of 
a large number of colors N c . In finite-temperature QCD, 
the fact that the number of gluon degrees of freedom is 
O(N^) while that of hadrons is 0(N®) hints at a first or- 
der confinement-deconfinement phase transition. At very 
large \ii thermodynamic quantities such as the density of 
isospin nj are proportional to N c . On the other hand, 
in the large N c limit the pion decay constant scales as 

= 0(N C ), and according to Eq. (0) the isospin density 
in the pion gas is also proportional to N c . Physically, the 
repulsion between pions becomes weaker as one goes to 
large N c , thus more pions are stacked at a given chem- 
ical potential. As a result, the N c dependence of ther- 
modynamic quantities is the same in the quark and the 
hadronic regimes. 

Confinement. — At large /zj, gluons softer than A are 
not screened by the Meissner or by the Debye effect [^5) : 
the condensate does not break gauge symmetry (in con- 
trast to the color superconducting condensate J3|) and 
there are no low-lying color excitations to screen the elec- 
tric field. Thus, the gluon sector below the A scale is de- 
scribed by pure gluodynamics, which is confining. This 
means there are no quark excitations above the ground 
state: all particles and holes must be confined in colorless 
objects, mesons and baryons, just like in vacuum QCD. 
If there is no transition along the /ij axis, we expect con- 
finement at all values of pLj. Since the running strong cou- 
pling a s at the scale of A is small, the confinement scale 
Aq CD (which is, in general, different from Aqcd) is much 
less than A. At large /ij, we thus predict a temperature 
driven deconfinement phase transition at a temperature 
T' c of order Aq CD , which is expected to be of first order 
as in pure gluodynamics. Since Aq CD <C A the hadronic 
spectrum is similar to that of a heavy quarkonium, with 
A playing the role of the heavy quark mass. 

The (T,fij) phase diagram. — By considering nonzero 
fij, we make the phase diagram of QCD three dimen- 
sional: (T, fiB, f^i). Two planes in this three-dimensonal 
space are of a special interest: the hb = (T,/ij) plane, 
which is completely accessible by present lattice tech- 
niques, and the T = (/i/,^B) plane, where the neu- 
tron star matter belongs. Two phenomena determine 
the (T, /i/) phase plane (Fig .[!]): pion condensation and 
confinement. 



3 



T 



:uy 5 d>=0 



<n->^0 <uy 5 d>*0 

m n |H/| 
FIG. 1. Phase diagram of QCD at finite isospin density. 

At sufficiently high temperature the condensate (||) 
melts (solid line in Fig. [l]). For large ///, this critical 
temperature is proportional to the BCS gap ([To|). There 
are two phases which differ by symmetry: the high tem- 
perature phase, where the explicit flavor U{1)l+r sym- 
metry is restored, and the low-temperature phase, where 
this symmetry is spontaneously broken. The phase tran- 
sition is in the 0(2) universality class Jig . The critical 
temperature T c vanishes at /i/ = and is an increasing 
function of /ij in both regimes we studied: |/xj| -C m p 
and \fij\ ^> Aqcd- Thus, it is likely that T c {fii) is a 
monotonous function of /i/. In addition, at large /x/, 
there is a first order deconfincmcnt phase transition at 
T' c much lower than T c (fii). Since there is no phase tran- 
sition at fii — (for small m u ,d) or at T = (assuming 
quark- hadron continuity), this first-order line must end 
at some point A on the (T,fj,j) plane (Fig. [I]). The exact 
location of A should be determined by lattice calcula- 
tions; one of the possibilities is drawn in Fig. |l]. 

The (/i/,/is) phase diagram. — This phase diagram de- 
serves a separate study. Here we shall only consider the 
regime ^> /is (the opposite limit /is 3> was con- 
sidered in Rcf. Jl?]]). Finite /ib provides a mismatch be- 
tween u and d Fermi spheres, which makes the supercon- 
ducting state unfavorable at some value of [Ib of order A. 
It is known Jr| that the destruction of this state occurs 
through two phase transitions: one at \±b slightly below 
A/y/2 and another at [j>b = 0.754A. The ground state 
between the two phase transitions is the Fulde-Ferrell- 
Larkin-Ovchinnikov (FFLO) state [fl8|| , characterized by 
a spatially modulated superfluid order parameter {v/y^d) 
with a wavenumber of order 2^. The FFLO state has 
the same symmetries as the inhomogeneous pion con- 
densation state which might form in electrically neutral 
nuclear matter at high densities It is conceivable 

that the two phases are actually one, i.e., continuously 
connected on the (///, /ib) phase diagram. 

The authors thank L. McLerran, J. Kogut, R. Pisarski, 
and E. Shuryak for discussions, the DOE Institute for 
Nuclear Theory at the University of Washington for its 
hospitality, and K. Rajagopal for drawing their attention 
to Ref. lil. 



[2] 
[3] 



I. Barbour et al, Nucl. Phys. B275 [FS17], 296 (1986); 
J.B. Kogut, M. P. Lombardo and D.K. Sinclair, Phys. 
Rev. D 51, 1282 (1995); Nucl. Phys. B, Proc. Suppl. 42, 
514 (1995). 

M.A. Stephanov, Phys. Rev. Lett. 76, 4472 (1996). 
B.C. Barrois, PhD Thesis, California Institute of Tech- 
nology, 1979; D. Bailin and A. Love, Phys. Rep. 107, 325 
(1984); M. Alford, K. Rajagopal and F. Wilczek, Phys. 
Lett. B422, 247 (1998), Nucl. Phys. B537, 443 (1999); 
R. Rapp, T. Schiifer, E.V. Shuryak and M. Velkovsky, 
Phys. Rev. Lett. 81, 53 (1998); Ann. Phys. (N.Y.) 280, 
35 (2000). 

M. Alford, A. Kapustin and F. Wilczek, Phys. Rev. D 
59, 054502 (1999). 

E. Dagotto, F. Karsch, and A. Moreo, Phys. Lett. 169B, 
421 (1986); E. Dagotto, A. Moreo, and U. Wolff, Phys. 
Rev. Lett. 57, 1292 (1986); Phys. Lett. B 186, 395 
(1987); S. Hands, J.B. Kogut, M.-P. Lombardo, S.E. 
Morrison, Nucl . Phys. B558, 3 2 7 (1999); S. Han ds and 



[6] 
[7] 
[8] 

[9] 

[10] 
[11] 

[12] 
[13] 
[14] 
[15] 

[16] 



[17] 
[18] 



[19] 



S.E. Morrison, |hep-lat/9902012j ; |hep-lat/9905021 
J.B. Kogut, M.A. Stephanov, and D. Toublan, Phys. 
Lett. B 464, 183 (1999); J.B. Kogut, M.A. Stephanov, 
D. Toublan, J.J. Verbaarschot and A. Zhitnitsky, Nucl. 
Phys. B582, 477 (2000). 

D. Weingarten, Phys. Rev. Lett. 51, 1830 (1983); E. Wit- 
ten, Phys. Rev. Lett. 51, 2351 (1983); S. Nussinov, Phys. 
Rev. Lett. 52, 966 (1984); D. Espriu, M. Gross and J.F. 
Wheater, Phys. Lett. 146B, 67 (1984). 
It is important, as is the case for 1 = 1, that there is no 
disconnected piece after ip integration in (|^). The proof 
does not apply, for example, to a or 77 meson correlators, 

r = i,7s. 

See, e.g., H. Georgi, Weak Interaction and Modern Par- 
ticle Theory (Benjamin-Cummings, Menlo Park, 1984). 
D.T. Son, Phys. Rev. D 59, 094019 (1999). 
WE. Brown, J.T. Liu, and H.-C. Ren, Phys. Rev. D 61, 
114012 (2000); 62, 054016 (2000); 62, 054013 (2000). 
A.J. Leggett, J. de Phys. 41, C7-19 (1980); P. Nozieres 
and S. Schmitt-Rink, J Low Temp. Phys. 5 9, 195 (1985). 
See, e.g., M. Randeria, cond-mat/9710223 and references 
therein. 

T. Schafer and F. Wilczek, Phys. Rev. Lett. 82, 3956 
(1999). 

This is similar to the behavior of the unbroken 
SU(2) C sector of two-flavor color superconductors (see 
D.H. Rischke, Phys. Rev. D 62, 034007 (2000)). 
The width of the Ginzburg region is suppressed by 
(A/ni) 4 at large fii and also by at large N c : 

J.B. Kogut, M.A. Stephanov and C.G. Strouthos, Phys. 
Rev. D 58, 096001 (1998) 



M. Alford, J. Bowers and K. Rajagopal, |hep-ph/000820£ . 
P. Fulde and A. Ferrell, Phys. Rev. 135, A550 (1964); 
A.I. Larkin and Yu.N. Ovchinnikov, Sov. Phys. JETP 
20, 762 (1965); see also Ref. @. 

A.B. Migdal, Sov. Phys. JETP 36, 1052 (1973); 
R.F. Sawyer, Phys. Rev. Lett. 29, 382 (1972); 
D.J. Scalapino, Phys. Rev. Lett. 29, 386 (1972). 



4 



